3 031 Hi N95- 10568 


Characteristics of a Drv. Pulsating 
Microburst at Denver Stapleton Airport. 


F. Proctor, 

NASA Langley Research Center 

[For these Proceedings, the author has furnished three papers: 

Influence of Low-Level Environmental Shear on Microburst Structure: 

Numerical Case Study, 

Numerical Simulation of a Pulsating, Low-Reflectivity Microburst Event, 

and 

Case Study of a Low-Reflectivity Pulsating Microburst: Numerical Simulation of 

the Denver, 8 July 1989, Storm. 

These appear in order on the following pages.] 
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INFLUENCE OF LOW-LEVEL ENVIRONMENTAL SHEAR 
ON MICROBURST STRUCTURE: NUMERICAL 
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Introduction 


In previous studies using the two-dimensional axisymmetric version of the Terminal 
Area Simulation System (TASS) the effect of environmental parameters and other factors 
was examined for isolated and stationary microbursts (Proctor 1988, 1989). The intensity 
and structure of the simulated microbursts were strongly affected by the lapse rates of 
ambient temperature and humidity, as well as the height of the melting level. These 
experiments, however, were limited to an axisymmetric model framework, and could not 
allow for an examination of the effects of translation and ambient wind shear. 

The present study extends these studies by examining the influence of ambient 
vertical wind shear on microburst intensity, asymmetry, and translation. This is achieved 
with the three-dimensional, high-resolution, version of the TASS model, by selectively 
varying the initial inputs for ambient low-level vertical wind shear and keeping all other 
input parameters fixed. A similar approach has been used by Weisman and Klemp (1982, 
1984) and Weisman (1993) in the study of vertical wind shear on isolated convective cells. 
Although our study here differs significantly from these studies in that we initially start out 
with a prescribed precipitation shaft as in Proctor (1988, 1989), rather than trying to 
simulate the complete storm evolution. 


Numerical Model 

Details of the numerical model used in this study have been described by Proctor 
(1987a), and should be referred to for details. The TASS model is a three-dimensional, 
time dependent, cloud model, with prognostic equations for momentum, pressure, 
continuity of water substance, and potential temperature. The six coupled equations for 
continuity of water substance include an equation each for water vapor, cloud droplets, 
cloud ice crystals, rain, hail or graupel, and snow. The model also includes 
parameterizations for surface friction and cloud-microphysical interactions. Salient 
characteristics of the TASS model are listed in table 1, and a list of the cloud 
microphysical parameterization are given in table 2. The TASS model has been validated 
against a number of cumulonimbus cases (e.g., Proctor 1987b, 1993; Proctor and Bowles 
1992), many of which were associated with intense microburst events. 


Initial Conditions and Description of Experiments 

An isolated microburst event is simulated by specifying a distribution of 
precipitation at the model top boundary and allowing it to fall within an environment 
conducive to wet microbursts. Twelve numerical experiments are conducted, which assume 
two precipitation shaft sizes and various ambient wind profiles as will be discussed below. 

The dimension of the physical grid domain is 10 km in the west-east direction (x- 
coordinate), 7.5 km in the south-north direction (y-coordinate), and 5 km in the vertical 
direction (z-coordinate). The domain is resolved by 32 levels with a vertical grid separation 
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stretching from 23 m near the ground to 280 m near 5 km. Each level is resolved by 102 
x 77 grid points with a uniform horizontal grid spacing of 100 m. 

The ambient (initial) profile for temperature and humidity typifies a warm moist 
wet-microburst environment, and is used for all experiments. 

The initial ambient winds are horizontally uniform and are a function of height only. 
The ambient winds are specified according to: 

U,(z) = V, * y c (^)°' 4 for z < 2000 m 

U e (z) = U 0 + U G for z * 2000 m 

and, 

Vfz) = V 0 

where U e and V e are the ambient wind components (m/s) in the x and y directions, 
respectively; and z is the altitude above the ground (m). At the ground (z = 0) the wind 
velocity is required to be zero in accordance with the no slip surface boundary condition. 
In all experiments the ambient wind speed and direction are constant above 2000 m. Input 
values for U G , U 0 , and V 0 are listed for each experiment in Table 3. 

For all experiments the type of precipitation prescribed at the top boundary is 
assumed to be of the hail/graupel class with a size distribution intercept of 2 x 10 5 m 4 . 
The distribution of precipitation at the top boundary is specified as 

Q(X,Y) = Q 0 e-w°? 


where 


R(X,Y) = JiX-XJ 2 + (Y-YJ 1 

and where X 0 ,Y 0 is a center location on the top boundary. For all experiments in this 
study Q 0 is 6 g m‘ 3 . The input parameter a, which controls the radial distribution of the 
precipitation shaft and ensuing downdraft diameter, is listed for each experiment in Table 
3. In all experiments the model grid and precipitation at the top boundary translate with 
the ambient wind velocity above 2000 m. 
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Table 1. Salient Characteristics of TASS 2.4 



Compressible, non hydrostatic equation set 
Non-Boussinesq formulation for density variations 
Three-dimensional staggered grid with stretched vertical spacing 
Movable, storm-centering mesh 

Explicit time-split, second-order, Adams-Bashforth time differencing and 
second-order quadratic-conservative space differencing for velocity and 
pressure 

Fourth-order quadratic-conservative space differencing and third-order 
Adams-Bashforth time differencing for temperature and water-vapor 

Third-order time/space differencing with upstream-biased quadratic 
interpolation for liquid and frozen water substance equations 

Nonreflective radiation boundary conditions applied to open lateral boundaries 

Filter and Sponge applied to top four rows in order to diminish gravity wave 
reflection at top boundary 

No explicit numerical filtering applied to interior points 
Surface friction layer based on Monin-Obukhov Similarity theory 
Smagorinsky subgrid-turbulence closure with Richardson number dependence 
Liquid and ice-phase microphysics 

Inverse-exponential size distributions assumed for rain, hail/graupel, and snow 
Raindrop intercept function of amount of rainwater 
Snow treated as spherical, low-density graupel-like snow particles 
Wet and dry hail growth 

Accumulated precipitation advected opposite of grid motion, so as to remain 
ground relative 

Radar reflectivity diagnosed from model rain, snow, and hail/graupel fields 
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Accretion of cloud droplets by rain 

Condensation of water vapor into cloud droplets 

Berry-Reinhardt formulation for autoconversion of cloud droplet 
water into rain 

Evaporation of rain and cloud droplets 

Spontaneous freezing of supercooled cloud droplets and rain 

Initiation of cloud ice crystals 

Ice crystal and snow growth due to riming 

Vapor deposition and sublimation of hail/graupel, snow, and cloud 
ice crystals 

Accretion by hail/graupel of cloud droplets, cloud ice crystals, rain, 
and snow 

Contact freezing of supercooled rain resulting from collisions with 
cloud ice crystals or snow 

Production of hail/graupel from snow riming 

Melting of cloud ice crystals, snow, and hail/graupel 

Shedding of unfrozen water during hail wet growth 

Shedding of water from melting hail/graupel and snow 

Conversion of cloud ice crystals into snow 

Accretion by snow of cloud droplets, cloud ice crystals, and rain 

Evaporation or vapor condensation on melting hail/graupel and snow 


183 




TABLE 3. Input variables for each experiment || 

RUN # 

u G 

(m/s) 

u 0 

(m/s) 

V 0 

(m/s) 

(J 

(m) 

AU e 

(m/s) 

FILE 

1 

0.0 

0 

0 

750 

0.0 

05a 

2 

5.0 

0 

0 

750 

4.2 

09a 

3 

12.5 

0 

0 

750 

10.4 

08a 

4 

17.5 

0 

0 

750 

14.6 

10a 

5 

25.0 

0 

0 

750 

20.8 

06a 

6 

0.0 

17.5 

0 

750 

0.0 

20a 

7 

17.5 

0 

8.75 

750 

14.6 

19a 

8 

0.0 

0 

0 

1350 

0.0 

15a 

9 

5.0 

0 

0 

1350 

4.2 

13a 

10 

12.5 

0 

0 

1350 

10.4 

12a 

11 

17.5 

0 

0 

1350 

14.6 

1 la 

12 

25.0 

0 

0 

1350 

20.8 

14a 


Values for the input parameters of the 12 experiments are listed in Table 3. As 
indicated by the value assumed for a, experiments 1-7 represent the small-diameter 
downdraft cases, while experiments 8-12 represent the large-diameter cases. The parameter 
AU e represents the vertical change in ambient wind from the lowest model level (at z = 23 
m) to the gradient level located at 2000 m, and thus is an indicator of the magnitude of 
ambient vertical shear. The ambient shear vector is directed eastward in all experiments 
(except Run #: 1, 6, and 8 which have no shear), and vanes in magnitude between each 
experiment. However, as in Run 7, the ambient wind vector may not always be in the 
same direction as the shear vector. The sensitivity of the microburst to ambient wind shear 
is examined by comparing the results between each of the simulations. 
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Results 


Sensitivity of key microburst parameters to environmental shear are shown in the 
following three line plots. The tendencies of both the small and large diameter microburst 
were similar and thus only results from the small diameter microburst are shown in the first 
two plots. 

The plot of ambient vertical shear vs peak velocity differential shows that increasing 
ambient shear reduces the peak velocity differential of a microburst. More importantly 
though, the peak velocity differential that is orthogonal to the ambient shear vector (e.g. 
AV along north-south segments) are less affected by increasing shear than the velocity 
differential along the shear vector (e.g. AV along east-west segments). Asymmetry of the 
microburst outflow field is indicated the by the difference in magnitude between the E-W 
and N-S segments. [For an axisymmetric outflow the peak E-W AV is equal to the peak 
N-S AV.] The results show that asymmetry of the outflow velocity field increases with 
increasing ambient shear. 

The second line plot shows the sensitivity of peak Fbar (1-km averaged F-Factor) 
at 280 m AGL. As described in Proctor and Bowles (1992) and Switzer et al (1993), 
values for F-factor are computed assuming horizontal north-south (or east-west) trajectories 
with an assumed air speed of 75 m/s. Similar to the peak AV curves in the previous plot, 
peak Fbar along segments orthogonal to the shear vector are stronger than the peak values 
along segments parallel to the shear vector, with widening differences as the magnitude of 
shear increases. But unlike in the previous plot, an optimal value occurs for Fbar along 
segments orthogonal to the shear vector. This optimal value happens for the north-south 
Fbar when the magnitude of the ambient shear is approximately 15 m/s. 

The third plot show the sensitivity of peak outflow speed (ground relative) for both 
large and small diameter microbursts. The peak outflow speed indicates the severity of 
the microburst in terms of its potential for inducing damage due to high winds. Our set 
of experiments indicates that the highest wind speeds occur for an optimal ambient vertical 
shear of about 15 m/s for both large and small diameter microbursts. The simulations 
indicate that this peak horizontal wind speed occurs on the downshear side of the 
microbursts. 

The last series of plots show cross sections from a case with strong ambient shear. 
The plots show radar reflectivity and velocity fields from Run 7 at a time near peak 
intensity. The ambient shear vector is directed east, although the storm is translating 
toward the east-northeast (see table 3). The figures show that the outflow is elongated E-W 
along the shear vector. The horizontal cross section of radar reflectivity at low levels 
shows a bow-echo shape with maximum radar reflectivity located on the downshear side 
of the outflow. The east-west vertical cross section shows a strong vortex circulation on 
the downshear side, which is associated with strong surface-level winds. Such a microburst 
as in this simulation, may be referred to as what Fujita (1985) calls a "rotor microburst." 
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SUMMARY AND DISCUSSION 


These results show that microburst asymmetry is influenced by the magnitude of 
the low-level ambient vertical shear. The microburst outflow elongates in the direction of 
the shear vector (which is not necessarily in the direction of translation), and generates the 
greatest hazard (for commercial jet transpons) along paths orthogonal to the shear vector. 
The model results also show that the asymmetry increases with increasing shear magnitude. 
One implication of these results concerns the detection of a microburst by a ground-based 
doppler systems. These systems may underestimate the hazard for landing and departing 
aircraft that are on trajectories orthogonal to both the sensor beam and shear vector, 
especially if the magnitude of the shear is large. Another implication is that microburst are 
more likely to be asymmetrical in regions (seasons) where there is climatologically a 
significant low-level shear. 

The model results also show that rotor microbursts and severe wind damage can be 
a product of the microburst interaction with strong ambient wind shear. 
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TERMINAL AREA SIMULATION SYSTEM 

(TASS) 

[ALSO KNOWN AS THE NASA WINDSHEAK MODEL] 


o 3-D TIME DEPENDENT EQUATIONS FOR 
COMPRESSIBLE NONHYDROSTATIC FLUIDS 

o PROGNOSTIC EQUATIONS FOR 11 VARIABLES 

1. 3-COMPONENTS OF VELOCITY 

2. PRESSURE 

3. POTENTIAL TEMPERATURE 

4. WATER VAPOR 

5. LIQUID CLOUD DROPLETS 

6. CLOUD ICE CRYSTALS 

7. RAIN 

8. SNOW 

9. HAIL/GRAUPEL 

10. INSECTS 

o lst-ORDER SUBGRID TURBULENCE CLOSURE 
WITH RICHARDSON NUMBER DEPENDENCY 

o SURFACE FRICTION LAYER BASED ON MONIN- 
OBUKHOV SIMILARITY THEORY 

o Choice At LATERAL BOUNDARIES Between OPEN 
BOUNDARY CONDITIONS ALLOWING MINIMAL 
REFLECTION or PERIODIC CONDITIONS 

o BULK PARAMETERIZATIONS OF CLOUD 
MICROPHYSICS 
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TERMINAL AREA SIMULATION SYSTEM 


CAN SIMULATE CONVECTIVE STORMS WITH REASONABLE 
COMPARISON TO REAL-WORLD EVENTS 

O REQUIRES INPUT VERTICAL PROFILE OF 
TEMPERATURE, HUMIDITY, WIND SPEED AND 
DIRECTION, AND ALTITUDE (OR PRESSURE) 

O SOUNDING MUST BE REPRESENTATIVE OF 
STORM’S ENVIRONMENT 

O SOUNDING MAY BE EITHER OBSERVED, 
FORECASTED ( e.g . MASS), OR INTERPOLATED 


SIMULATIONS FROM TASS EXTENSIVELY VALIDATED 
AGAINST OBSERVED STORMS RANGING FROM 

O SHORT-LIVED WEAK CUMULONIMBUS 
O SEVERE LONG-LASTING SUPERCELL STORMS 
O MICROBURST PRODUCING STORMS 
O HAILSTORMS AND TORNADIC THUNDERSTORMS 
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MICROBURST SENSITIVITY TO AMBIENT VERTICAL SHEAR 
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WET MICROBURST ENVIRONMENT 
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MICROBUST PEAK WIND DIFFERENTIAL AND PEAK F-FACTOR 
ARE FOUND IN DIRECTION NORMAL TO SHEAR VECTOR 



Asymmetric Microburst - Run 7 

Environmental Shear Directed Toward East 
Storm Directed Toward ENE 
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Wind Speed at 50 m AGL 
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X (km) Zfoyf 

Radar Reflectivity and Wind Vectors at 50 m AGL 


Asymmetric Microburst - Run 7 

Vertical Cross-Section of Radar Reflectivity and Wind Vectors 
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Storm Relative Wind Vectors 
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NUMERICAL SIMULATION OF A PULSATING, 
LOW-REFLECTIVITY MICROBURST EVENT 
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SUMMARY 
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Denver, 8 July 1989, Simulation 
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DEVELOPMENT OF THE STORM 
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Denver, 8 July 1989, Microburst Simulation 

North-South F-Factor at 44 min and 1 50 m AGL 
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Denver, 8 July 1989, Microburst Simulation 
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SOME LOW-REFLECTIVITY MICROBURST EVENTS MAY BE DIFFICULT 
TO DETECT WITH RADAR DUE TO LACK OF BOTH PRECIPITATION 
AND INSECTS 



INSECT MODEL 
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ave a maximum fall speed of 75 cm/s 



Local Number Concentration of Insects Governed by: 
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Based on insect behavior as reported in: G. L. Achtemeier; The Use of Insects as Tracers for "Clear- 
Air" Boundary-Layer Studies by Doppler Radar, J. Atmos Oceanic Technol.. 8, 1991, pp 746-765. 



